Engineering antiferromagnetic Heisenberg spin chains for maximizing of the 

groundstate entanglement 
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We study the correlation function and concurrence for the eigenstates with zero spin of engineered 
Heisenberg models to explore the entanglement property. It is shown that the total nearest neighbor 
(NN) correlation function of zero-spin eigenstates reaches its local extremum when the coupling 
strength is uniform, and correspondingly the groundstate entanglement of d-T) cubic AF Heisenberg 
model is locally maximized. Moreover, numerical calculations for a N-siie quantum spin ring with 
cosinusoidally modulated exchange coupling, i.e. Ji = J(l + cos(27rz/iV)), indicate that the uniform 
coupUng is not the unique optimal distribution for maximizing the groundstate entanglement and 
this modulation of interactions can induce the longer range entanglement. 
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Introduction. Studies of various quantum spin mod- 
els are very important in understanding the properties 
of realistic solid state systems. Numerical and analytical 
studies have shown that the groundstate properties are 
sensitive to the geometry of the lattice 111 |^ y, |5{ 
and the distribution of the exchange couplings between 
two neighbor spins 0, S IS 0- Although much effort 
has already been made to investigate the properties of 
the ground state, the exact result about it is still rare. 
For instance, a spin-1/2 Heisenberg ring is one of the 
nontrivial, relatively simple physical systems. It was the 
first one dimensional quantum model solved exactly by a 
straightforward diagonalization of the Hamiltonian and 
has been used as a testing ground for many theoretical 
approaches. Bet he found the eigenvalues and eigenfunc- 
tions of the corresponding Hamiltonian 0, but the ob- 
tained exact eigenfunctions are so complicated that they 
are almost useless for calculating some physical quanti- 
ties, such as something relevant to correlation functions, 
to reveal all the characteristic properties of the ground 
state. Recently, quantum entanglement was proposed as 
a new type of physical resource, which is desired to de- 
pict the features of the ground state, like the correlation 
function or other conservative quantities. It was found 
that the ground state of an AF Heisenberg ring possesses 
maximal entanglement compared to that of excited states 
[m and the entanglement is believed to have somthing 
to do with the quantum phase transitions [l^. 

The present paper explores the pairwise entanglement 
of the zero-spin ground state of the Heisenberg model 
with certain engineered inhomogeneous coupling con- 
stants. We first make a simple observation that the total 
nearest neighbor (NN) correlation function of zero-spin 
eigenstate (with vanishing components of total spin) is 
locally maximized for the uniformly distributed coupling 
strengths, while the groundstate entanglement of the d- 
D cubic antiferromagnetic (AF) Heisenberg model can 
locally maximize in this case, correspondingly. How- 
ever, the uniform distribution of coupling constants is 



not the unique optimal one for maximizing the ground- 
state entanglement in the case with engineered coupling 
constants. Namely, the homogeneity of couplings is only 
the sufficient conditions rather than a necessary condi- 
tion in maximizing the groundstate entanglement. To 
demonstrate this point of view, we study the A/'-site ring 
system with cosinusoidally modulated exchange coupling, 
i.e., Ji = J(l + cos{27Ti / N)) in both analytical and nu- 
merical approaches. We discovered that the ground state 
with varying couplings Ji and that with the fixed cou- 
pling J are nearly identical. 

Maximization of the groundstate entanglement. The 
Hamiltonian H = ^{^ij^ Jij^i • of Heisenberg spin 
model can be written in terms of the spin operator at 
ith site, where Jij is the coupling constant of exchange 
interaction, which is trivially restricted to be nonzero in 
this paper. Our studies will focus on those states with 
vanishing components of total spin, or called the zero 
spin states. In this sense the relationship between corre- 
lation and concurrence has been well established 13]. To 
describe the behavior of correlation one can define the 
nearest neighbor (NN) correlation function 



(1) 



(ij) 



for the eigenstate |?/^o) with zero spin. The sandwich 
{Si • Sj)o = {^|Jo\Si ■ Sj|?/^o) is defined as the expectation 
value of Si ■ Sj in the zero-spin eigenstate \ipo). For- 
mally, it is proportional to the average of the Hamiltonian 
by assuming the identical coupling constants. What we 
concern is the J^j -dependent behavior of the correlation 
function, which can be characterized by the extremum of 
Fo{Jij). 

Now, for a zero-spin eigenstate, we show that Fo{Jij) 
can reach its extremum when all the exchange constants 
are identical, Jij = J, if Fo{Jij) is analytical at this point. 
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Actually, the eigen energy can be written as 



(2) 



Differentiating the above equation with respect to an ar- 
bitrary coupling strength Jki^ one can get dEo/dJki = 
S(ii> J^j^^^i ' ^j)o/dJki + (S/c • S^)o. On the other hand, 
together with the Feynman-Hellman theorem 



dEo _ dH 



(3) 



it indicates that Xl^^j) Jijd{Si ■ Sj)o/dJik = 0. Obvi- 
ously, when all the coupling strengths Jij = J, we get 
dFo/dJij = 0, which means that Fq has an extremum 
for the uniform coupling strength. Notice that this con- 
clusion is always true no matter the system is bipartite 
or non-bipartite lattice and the state is ground state or 
excited state. However, Jij = J may not be the unique 
distribution for Fq to reach the extremum. We will find 
that there may exist a periodic coupling strength dis- 
tribution in 1-D system which induces the same ground 
state approximately as that with uniform coupling dis- 
tribution. 

Now we apply the above conclusion to the bipartite 
lattice with Jij > and Na = Nb^ where Na^ Nb are 
the numbers of the sites belongin g to sublattices A and 
B. According to Lieb's theorem 'l4], the ground state 
is singlet and has S = 0, which ensures that Fg reaches 
its minimum Fg{l) = Eg/N at the point Jij = 1, where 
the zero-spin state |?/^o) is replaced by the ground state 
1^). Furthermore, we consider the groundstate entangle- 
ment for the AF Heisenberg model on a d-D cubic lat- 
tice with translational symmetry in all directions. Here 
the symmetry specifies the geometry of the lattice only, 
i.e., the Hamiltonian may not have the translational sym- 
metry. Since the upper bound of Eg is — X^^^^-^ 1/4 in 
the vicinity of the point Jij = 1, we have the inequality 
Fg{l) < -(l/iV)E(i,) 1/4 from (S^ • S,)g < -1/4. On 
the other hand, the pairwise concurrence for such system 
is Cg = (l/2)max{-(4/iV)E(i,)({Si ■ S,), - 1),0} = 
(l/2)max{-[4F<,(l) + (1/iV) Efe) 1],0} M, or 

C, = -i[4F,(l) + ^], (4) 

which maximizes in the vicinity of the point Jij = 1, 
where ^ = ^ ^^^^ numbers. Therefore we get 

the conclusion that the ground state has locally maximal 
pairwise entanglement when the exchange interactions 
distribute uniformly. A similar conclusion has been ob- 
tained for the XX Z model at the isotropic point [l5L[l^ . 

Notice that the above statement does not mean that 
the concurrence is the maximum in the whole range of 
coupling constants, but just in the vicinity of the uniform 
point. In order to illustrate this, we investigate a sim- 
ple spin model, the AF Heisenberg ring with alternative 



coupling constant. The Hamiltonian of N-site ring is 



i^odd 



i^even 



It is well known that the value of the NN pairwise 
concurrence is 0.386 for J = 1 13]. On the other 
hand, in the limit cases of J = and oo, the exact 

ground states are (/)i = [12] [34] [N - IN], and 

(j)2 = [23] [45] [A/'l] respectively, where [ij] denotes 

a resonant valence bond (RVB) of two spins located at 
the lattice sites i and j 17]. Both and 02 has a same 
concurrence 0.5, which is larger than that of the state at 
uniform point. This demonstrates that the concurrence 
is not the maximum in the whole range of coupling con- 
stants. The groundstate concurrences for small systems 
with N = 8, 10 and 12 are calculated by exact diagonal- 
ization method and plotted in Fig.l. It indicates that 
the concurrences takes its local maximum at the uniform 
point. 

C sinus oidally modulated systems. As mentioned in 
advance, the uniform system may not be the unique op- 
timal one possessing a maximal groundstate entangle- 
ment. A simplest case is that there may exist a model 
with non-uniform coupling strength distribution, which 
has the same ground state as that of the uniform one. We 
now show this observation in the XY spin ring exactly 
and the XXX spin ring approximately. 
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FIG. 1: The groundstate concurrences vs the coupling strength 
for the alternative Heisenberg ring systems with N = 8 (cir- 
cle), 10 (solid circle) and 12 (solid line). It shows that the 
uniform coupling point is the local maximum. 

We consider the two cases of spin- 1/2 Heisenberg 
model with the anisotropy parameter A = 0, 1 on an A^- 
site ring. With the cosinusoidally modulated exchange 
couplings, i.e. = J(l + cos(27ri/A^)), the Hamiltonian 
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H = Hq -\- Hadd can be separated into two parts 

N 

Ho = J Y^iSfSf^, + SfSf^, + A5f (6) 



and 

TV . 

Hadd = J^cos(^)(5f5,^i + A^f^f^,). 

i 

(7) 

where 27Ti/N is a factor to determine the profile of the 
additional Hamiltonian. H was a pure toy model before 
the array of quantum dots is considered as a media to 
transfer the c Q Q Q " Fig.2 H and 

Hq are illustr see that Hadd 

has a subtle r nd states. 



(a) 



(b) 



FIG. 2: Schematic illustration of co sinus oidally modulated (a) 
and uniform (h) NN coupling spin ring systems. The two most 
separated sites in (a) is the nearest neighbor in (h). 

The simplest case A = is first considered as the 
so called XY model or the hardcore boson model. It is 
well known that if N/2 is odd the spectrum of the XY 
Heisenberg spin ring in the subspace with = is re- 
duced to the model of non-interacting spinless fermion 
system. The reduced Hamiltonian consists of two parts 
Hr = {JmY'^iala^+.+h.c) and H^,\ = (J/2) 
Y^f cos{27ii/N)(a\aij^i + /i.c), written respectively in 
terms of the spinless fermion operata a\ at zth site. In k- 
space, by using the discrete Fourier transformation a\ = 
alexp{ikl),{ k = 2iin/N,n = 0,±1,±2,..., 
±(A^/2 — 1), N/2)^ they can be re- written as 



H, 



XY 



J COS ka\ak 



and 



H, 



XY 
add 



(8) 



The physics of the above Hamiltonian is quite ob- 
vious. The additional Hamiltonian H-^^ describes an 
A^-site chain system with nearest neighbor hopping in 
k space. Notice that the /c-dependent hopping inte- 
gral (X cos(/c + j^) vanishes at points k = -^kf^ where 



kf=7T/2 — ir/Nis the "fermi point" ( surface ) for the 
one dimensional half- filled spinless fermion model H^^ . 
It is easy to find that the particles are confined either in 
the regions \k\ < kf or \k\ > kf. Actually, it is easy to 
prove that the particle number in the region \k\ < kf ( 
or \k\ > kf ) is a good quantum number for H^J^, i.e., 

[ n, i^^I ] = 0, where n = Eke{\k\<kf} 4^^. 

Obviously, there are only four eigenstates of Hq^ with 
zero eigenvalues: the groundstate li^f^) = n|/c|</c^ ^II^) 

and the eigenstate {^p^^J = Yl\k\>kf of ^o^ with 
maximum eigenvalue. These two states indicate that 
the two regions in /c-space separated by fermi points are 
fully filled respectively. If all the /c-space is fully filled 
or empty, one can get other two eigenstates |V^fm), be- 
ing the saturated ferromagnetic states with all spin up 
or down. These four eigenstates are also the eigenstates 
of any linear combinations of Hq^^ and H^^. Here we 
only consider the simplest case H^^ 
that is 



J^XY I rrXY 
^0 ^ ^add ' 



H^^ = 2jf2cosH^){ala,+, + h.c). 



(9) 



When N/2 is even, the XY model is also equivalent to 
the hardcore boson system. But we can not get the same 
analytical result as that in the case of odd number of 
N/2. Although numerical results for small N-site sys- 
tem show that the same conclusion is also true, we can 
not give an exact proof at present stage. Anyway, this 
sample implies that the factor 27Ti/N is a characteristic 
groundstate property for the spin ring system. 

Finally we turn to the isotropic case, i.e. A = 1, in 
which the Hamiltonians read as Hq = Si • S^+i 

and Hadd = J^f cos{27ri/N)Si ■ S^+i. Correspondingly 
we have 



N 



H = Ho^ Hadd = 2J^cos2(^)S, • S,+i. (10) 



For J < 0, it is easy to find that the ground states of Hq 
and H are both saturated ferromagnetic with the same 
eigen energy NJ /A. The saturated ferromagnetic state is 
also the eigenstate of Hadd with zero eigenvalue. 

For J > 0, the saturated ferromagnetic state is also 
the common eigenstate of i^o, Hadd and H. So part 
of the exact conclusion obtained with XY model can 
be extended to the isotropic antiferromagnetic Heisen- 
berg spin ring trivially. For the ground state, we can 
not get similar analytical result as that in XY model. 
In order to investigate the relation between the ground 
states of H and i^o, exact diagonalization is performed 
to compute the groundstate energies and the overlap of 
the two groundstate wavefunctions. In Table 1, the nu- 
merical results are listed for the systems of site number 
N = 12, 14, 16, 18, 20, 22 and 24. It shows that the over- 
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lap approaches the unity, i.e. the ground states of these 
two Hamiltonians are approximately identical. 



N 






Overlap 


12 


2.5x10" 


-5 


0.999988 


14 


2.6x10" 


-5 


0.999983 


16 


2.5x10" 


-5 


0.999978 


18 


2.4x10" 


-5 


0.999974 


20 


2.3x10" 


-5 


0.999971 


22 


2.2x10" 


-5 


0.999967 


24 


2.0x10" 


-5 


0.999964 



Table 1 

Table 1. The differences of groundstate energies 
and the overlaps of corresponding eigenf unctions of the 
Hamiltonian H and Hq for the N-site systems obtained 
by exact diagonalization. It shows that the ground states 
of the two Hamiltonians are approximately identical. 
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FIG. 3: The overlaps between the groundstate wavef unctions 
of Ho and Hadd with the coupling strength Ji = J(l + 
0.5cos{2n7ri/N))(a) and Ji = J(l + 0.95 cos(2n7ri/A^)) (b) 
for N = 20^ 18 and 16. It shows that only 27Ti/N matches the 
ground state of uniform coupling system. 

The above results imply that the homogeneity of cou- 
plings is not the unique optimal distribution for the 
translational invariant groundstate having maximal en- 
tanglement. Now we investigate the ground states of 
the N-site ring systems with the coupling strength Ji = 
J ^ J' cos(2n7ri/7V)), where n = 0, 1, 2, N - I. Here 
J' < J ensures > in order to avoid the degeneracy 
of the ground states 14]. Obviously, for XY model the 
conclusion for ground state obtained above is no longer 
available for n ^ 0, 1. For isotropic Heisenberg model, a 
small size ring is investigated by numerical method. In 
Fig. 3, the overlap of the corresponding eigenfunctions are 
plotted. It shows that 2i\ jN is the characteristic factor 



for the isotropic Heisenberg ring system. On the other 
hand, from Fig. 2 we can see that the concurrence of spins 
at the two ends of a modulated chain should be close to 
that of NN spins on a uniform ring. 

Summary. In summary, the correlation and entan- 
glement of the zero-spin eigenstates of the Heisenberg 
models are studied. We find that the total NN corre- 
lation function of zero-spin eigenstate reaches its local 
extremum when all coupling strengths are identical. Ap- 
plying this fact to d-D cubic AF Heisenberg model, the 
groundstate concurrence, the measure of entanglement is 
locally maximized at the same point. Numerical calcula- 
tions are employed to investigate a N-siie quantum spin 
ring with cosinusoidally modulated exchange couplings. 
It indicates that the homogeneity of couplings is not the 
unique optimal distribution for maximizing the ground- 
state entanglement and this modulation of interactions 
can indeed result a longer rang entanglement. Further- 
more, it also implies that 27t/N is the characteristic factor 
for the isotropic Heisenberg ringe system. 
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